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Based on ARPES and STM work, Co doped NaFeAs has been reported to be a strongly coupled 
superconductor, similar to cuprates However, a recent ARPES result, claims otherwise, with 

2A/kBTc slightly beyond the BCS hmit, 3.53 [3|- We have performed ^^As and "Na nuclear 
magnetic resonance(NMR) on high-quality, single crystals of NaFeo.975Coo.025 As. The temperature 
dependence of the '^^As Knight shift, fitted to a Yosida function, suggests moderately weak-coupling 
superconductivity in a singlet state with 2A/fc_BT'c = 4.0 ± 0.3. The temperature dependence of 
the penetration depth obtained from ^^Na NMR was found to be consistent with the Knight shift 
results and Aa6(0) was found to be 5, 327 ± 78A. 

PACS numbers: 



Superconductors have been classified as weak or strong 
coupling depending on the amplitude of the energy gap 
relative to T^, where l^lkoTc = 3.53 is the weakcoupling 
BCS limit. Larger values can be associated with strong 
coupling and a possible indication of strong electronic 
correlations. '^He, cuprate superconductors, and heavy 
fermion materials are in this class. However, a consis- 
tent picture for the recently discovered pnictide super- 
conductors has not yet emerged with a wide spread in 
the reports of 2A/fcBTc Q from 1 to 8. K-subsitituted 
BaFe2As2 (hole doped) has doping independent behavior, 
2A//cbT'c ~ 3.1, found from point contact Adreev reflec- 
tion (PCAR), whereas Co-substituted BaFe2As2 (elec- 
tron doped) shows doping dependence with strong cou- 
pling superconductivity near optimal doping 
might be inferred, therefore, that Co-substituted NaFeAs 
(electron doped) should also be strongly coupled, and in- 
deed this was reported from angle resolved photoemis- 
sion spectroscopy (ARPES) and scanning tunneling mi- 
croscopy (STM) [l-d However, recent ARPES ex- 
periments question these claims arguing that the super- 
conducting gap size had been over-estimated owing to 
larger impurity scattering rates Q. Since ARPES and 
STM measurements are surface sensitive, it is impor- 
tant to investigate these issues with a probe that de- 
tects bulk superconductivity. NMR satisfies these re- 
quirements, and here we report that superconductivity 
in NaFeo.975Coo.o25As (NaCo25) has moderately weak- 
coupling from the temperature dependence of the Knight 
shift and penetration depth in high-quality NaCo25, 
taken from measurements of ^^As and ^^Na NMR re- 
spectively. 

Our ^^As and ^^Na NMR experiments were performed 
at Northwestern University from 4.2 K to room temper- 
ature with external magnetic field H = 6.4 T parallel 
to the c-axis. Because of the high reactivity of Na we 
used an environmentally sealed sample holder. Fig. 1 (a), 
made out of glass, for "^^As NMR, or Stycast 1266® for 



^^Na. The containers were filled with He in an oxygen 
free chamber to prevent sample degradation and to ex- 
change heat. Degradation can be identified over time 
from an increase in the NMR linewidth and consequently 
different crystals were used for ^^As and ^^Na NMR. In 
the work we report here, both exhibited very narrow 
NMR spectra at room temperature: 8 kHz for ^^As, and 3 
kHz for ^^Na. Crystals with dimensions ~ 3 x 2 x 0.3 mm^ 
were grown at the University of Tennessee and found to 
have Tc of 21 K from magnetization with H = 10 G par- 
allel to the a&-plane, Fig. 1 (b). Typically, spin-echo 
sequences (7r/2 - tt) were used to obtain spectra. Knight 
shift, and spin-lattice relaxation rate, l/Ti, for the cen- 
tral transition (-1/2 ol/2) with a 7r-pulse w 8 usee. 
The spin-spin relaxation time, T2, was measured with 
a Carr-Purcell-Meiboom-Gill (CPMG) sequence and Ti 
was measured with the usual saturation recovery method 
when it was shorter than 2 sec. For longer Ti the more 
efficient progressive saturation technique was used. 

The "^^As and "Na NMR spectra from the central 
transitions are shown in Fig. 1 (c) and (d). The nar- 
row linewidths of the spectra in the normal state, 10 
kHz C^^As) and 3 kHz ("Na) at 20 K, and their very 
weak dependence on temperature are indications of high 
crystal quality. In the normal state, cooling from high 
temperatures, both ^^As and -'^-'^Na spectra shift to lower 
frequency as has been observed in other pnictide super- 
conductors [l3€l- Below Tc = 18 K in 6.4 T, the abrupt 
decrease in the ^^As spectra indicates a spin-singlet tran- 
sition. Formation of the vortex lattice is evident in broad- 
ening of the spectra in the superconducting state. 

The measured NMR frequency shift, Aoj{T), relative 
to the Larmor frequency, lS.uj{T) /ujl = (i-^iT) — uji,)/ujl, 
can be expressed as, 

K{T) = Ks{T) + Korb + Kq+ 47r(l - D)M{T)/H, (1) 

generically referred to as the Knight shift, ^^AT and 
^^K for each nucleus, determined from the peak posi- 
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FIG. 1. (a) Schematic drawing of the sample holder; (b) mag- 
netization measurements; (c, d) NMR spectra of NaCo25 in 
6.4 T with H || c-axis. The narrow linewidths confirm high 
sample quality. As with other pnictide superconductors [iol - 
[T^ NaCo25 shows a continuous decrease on cooling from room 
temperature above Tc- Below Tc = 18 K (blue trace) the spec- 
tra shift noticeably to lower frequency with decreasing tem- 
perature indicating spin-singlet superconductivity in NaCo25. 
The increased linewidth on cooling below Tc is due to the for- 
mation of a vortex lattice. 



tions of the spectra m Fig. 2. The Knight shift has 
electron spin, Kg, orbital, Korb, quadrupolar, Kq, and 
Kdia terms, the latter from the superconducting diamag- 
netism, M{T) with demagnetization factor, D; Kort and 
Kq provide a temperature independent offset. In the 
normal state the temperature dependence comes from 
Ks{T) = Xs{T)Ahf, given by the spin susceptibility, 
Xs{T), and its corresponding hyperfine coupling, Ahf- 
For r > 100 K the different temperature dependence 
of '^^K and ^^K is due to a much smaller hyperfine 
field for Na as compared to As, which we estimate to 
be Ahf 1'^'^ Ahf ~ 19. For T < 100 K the compar- 
ison oi ^^K with ''^if indicates ^^Kg is negligible and 
"i^orh +" Kq = 0.0033%. The decrease of "if in the 
superconducting state is due entirely to Kdia{T) and pro- 
vides an internal measure of the local field that can be 
subtracted from ^^/v to give a precise determination of 
the spin shift, '^^Ks{T < Tc), Fig. 2(b). 

We assume a single isotropic superconducting energy 
gap in NaCo25. Even though the sizes of the super- 
conducting gaps from different ARPES measurements in 
NaFeo.95Coo.05As (NaCo50) are different, they are con- 
sistent in two ways: no anisotropy in the superconduct- 
ing gaps in NaCoSO, and the sizes of the superconduct- 
ing gaps on the electron and hole Fermi surfaces are al- 
most identical [H, 01 • Consequently, we fit the tempera- 
ture dependence of '^^Kg to an isotropic Yosida function, 
Y{T) [ill, Eq. 2, to determine A(r). Wc take the gap to 
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FIG. 2. Knight shifts for ^''As and ^^Na over a wide tempera- 
ture range, (a) The decrease of K in the normal state on cool- 
ing is characteristic of pnictide superconductors, (b) The low 
temperature behavior of "^^Ka, after subtraction of Kdia, is 
fitted to a Yosida function solid curve, with 2A(0)/fcsTc 
= 4.0 which we interpret as evidence for moderately weak- 
coupling superconductivity. 



be scaled [3] to the weak-coupling BCS superconducting 
gap, Ao (T), represented in Eq. 3 by a convenient analyt- 
ical form [l5| , and then we fit the measured Knight shift 
to determine the scale factor a and hence A(0). 
For T < Tc, 



Ks{T) 
Ks{Tc) 

A{T) = aAo(T) 



Y{T) = 1 - 2TTkBTY^ 



A(r)2 



(4 + A(T)2)3/2 



(2) 



Ao(r).Ao(0)tanh(||^^|(|-l)f 



(3) 
) (4) 



where e„ = 27r(n -I- 1/2)^^^ ^ and ^|o is 1.43. 
With this procedure we have found the superconduct- 
ing gap size to be 3.6 ± 0.3 meV, 2A/ ksTf. = 4.0, and 
a = 1.14 ± 0.1. The resultant fit is shown in Fig. 2.(b). 
Our result is slightly larger than the BCS limit but still 
supports the basic idea that superconductivity in NaCo25 
is weakly-coupled consistent with a recent ARPES mea- 
surements Q ■ The Knight shift is a thermodynamic mea- 
surement and at low temperatures is sensitive to the low 
lying excitations < 1 meV above the energy gap. If in- 
deed there are two gaps of comparable size we should 
interpret our results in terms of the smaller of these. Re- 
cent STM measurements [1] on the same material, and 
from the same source as ours, infer a larger energy gap of 
5.5 meV. It is interesting that the bosonic mode reported 
in that work is within 8% in energy to 2A inferred from 
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the Knight shift. The disparity in energy gap size be- gives the open green triangles in Fig. 3. Subtracting 
tween Knight shift and STM may be a reflection of the this spin-fluctuation component of l/-^^TiT leaves a sec- 
different sensitivity in the two experiments to gap struc- ond contribution represented by the dashed red curve, 
ture and energy scales. We speculate that the spin-fluctuations come from g w tt 

and that the remainder might be from, q ~ 0, not signif- 
1 [ ' ' ' ' ~j icantly affected by superconductivity. 




FIG. 3. The spin-lattice relaxation rate divided by temper- 
ature, l/TiT, from ^^As and "Na. The increase of 
(blue circles) on cooling in the normal state is due to spin 
fluctuations and drops for T < Tc with the onset of su- 
perconductivity. Despite the small hyperfine field, 
has a similar behavior near Tc (red circles). Scaling the rates 
by the square of the ratio of the hyperfine flelds and their 
respective gyromagnetic ratios gives the open green triangles 
for the expected contribution at the Na site. What remains 
is a second component of relaxation (red dashed curve) of 
unknown origin. 

In contrast to the Knight shift (oj w 0, q « 0) the spin- 
lattice relaxation is determined by the sum over wave 
vectors, g, of the imaginary part of the dynamic sus- 
ceptibility which is coupled to the nucleus through the 
square of the hyperfine field given by the form factor, 
F{q), Eq. 5, 



l/TiTcx72^i^(<7) 



Im[x±iq,uj)] 



(5) 



where 7 is the gyromagnetic ratio, F[q) the form factor, 
Im [x±] is the imaginary part of the dynamic suscepti- 
bility perpendicular to the field. 

We have found different behavior for l/^^TiT and 
l/^TiT , Fig. 3. In the normal state, if'^TiT increases 
with decreasing temperature approaching 7^ and can be 
understood in terms of spin fiuctuations E, 12, 17, 18l | 
expressed in terms of the Moriya theory [19[, Eq.5. The 
increase is suppressed with the onset of superconductiv- 
ity, so that l/'^^TiT drops dramatically below Tc, ap- 
parently strongly dependent on the density of normal 
quasiparticle excitations. For Na we find that l/^^TiT 
has a component that mimics this behavior and can be 
accounted for by the weak hyperfine coupling strength, 
^^Ahf- Scaling the rates by the square of the ratio of the 
hyperfine fields and their respective gyromagnetic ratios, 



FIG. 4. Temperature dependence of the linewidth, full- 
width at half-maximum (FWHM), of the '^'^As and "Na cen- 
tral transition, and the spin-spin relaxation time, T2. The 
narrow hnewidths at T = 20 K from ^^As (10 kHz) and "Na 
(3 kHz) and the weak temperature dependences indicate high 
crystal quality. Below Tc the linewidth increases in the inho- 
mogeneous vortex state and becomes constant below « 13 K 
for ^^As (12 K for "Na). 



We next turn our discussion to the diamagnetism and 
vortex supercurrents in the superconducting state. The 
linewidths of the central transitions of ''^As and ^^Na 
are shown in Fig. 4. Below T,. the linewidths from 
both samples broaden due to vortices and approach a 
constant at low temperature approximately at the min- 
imum in ^^T2, which we identify with the formation 
of a robust vortex solid. In our earlier work on Co 
doped BaFezAsz [l^ we have associated a minimum T2 
with the irreversibility temperature where vortices be- 
come pinned. Although ^^Na NMR is not as effective 
as ''^As NMR for probing the electronic excitations, it 
nonetheless provides a higher resolution tool for analy- 
sis of the spectrum in the superconducting state which 
is affected by vortices and from which we can determine 
the penetration depth, Aqj,. The comparisons between 
the superconducting state spectra and the normal state 
spectra in Fig. 5, indicate an increase in the high fre- 
quency portion of the spectrum as expected in a vor- 
tex solid. We have solved the Ginzburg-Landau (GL) 
equations in an ideal vortex lattice using Brandt's algo- 
rithm [2l[, where we set i?c2 = 36 T (23 • The resulting 
best fit distribution of local fields for ^-'^Na at T = 4.2 
K is presented in Fig. 5 (c). The fit was performed 
by convoluting the GL solution from the algorithm and 
the normal state [T = 20 K) spectrum, minimizing its 
difference from the experimental spectra using the pene- 
tration depth as an adjustable parameter. We have an- 
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FIG. 5. Spectra in the normal state (20 K, shifted in fre- 
quency for comparison) and the superconducting state (4.2 
K) for (a) As and (b)"Na. The hnewidth broade nine due 
to vortices is evident in the superconducting state, especiaUy 
clear for '^^Na. (c) The Ginzburg-Landau calculation that 
best fits the superconducting state ^^Na spectra at T = K 
from which the penetration depth, Aai,, was determined, (d) 
Penetration depth as a function of temperature where Aa6(0) 
is 5, 327 ± 78 A. Temperature dependence of the penetra- 
tion depth is shown in (d) compared with a /iSR experiment, 
Aab(O) = 4,260l [131. Fitting to a BCS formula we found 
the superconducting gap to be 3.0 ± 0.3 meV consistent with 
Knight shift data. 



alyzed Xab{T) using a BCS formula in the low temper- 
ature limit, A(r) = A(0)[1 + y^i^exp(-^)] 

For ^5 As we find Xab{0) to be 4, 780 A and from "Na, 
XabiO) = 5, 327 A, which can be compared with fiSK poj . 
Xab{0) = 4, 260l. The penetration depth from the "Na 



is more reliable since the temperature dependent spin 
shift is negligible. From the BCS fit we found the su- 
perconducting gap size 3.0 ± 0.4 meV, close to the value 
from the Knight shift, supporting the view that NaCo25 
superconductivity is close to weak-coupling. 

In summary, we have found evidence for moderately 
weak-coupling superconductivity from the Knight shift 
and the penetration depth in NaCo25. The supercon- 
ducting gap sizes, 3.6 ± 0.3 meV, from the Knight shift 
and 3.0 ± 0.4 meV from the penetration depth are con- 
sistent with 2A/kBTc = 4.0 ± 0.3 taken from the Knight 
shift results. Additionally, the spectra in the mixed state 
have been successfully represented in Ginzburg-Landau 
theory, giving a penetration depth of Aa6(0) = 5,327± 
781. 
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